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OO Abstract 

in 
in 

^ Using integrability techniques, we compute four-point functions of single trace gauge-invariant 

,__! operators in A/" = 4 SYM to leading order at weak coupling. Our results are valid for oper- 

^""1 ators of arbitrary size. In particular, we study the limit in which two of the four operators 

>■ are taken to be much smaller than the others. We show that in this limit our weak coupling 

K^ result matches with the strong coupling result in the Frolov-Tseytlin limit. 
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1 Introduction 

Recently, there has been a lot of activity in computing three-point functions of single trace 
operators in the context of the AdS/CFT correspondence [1, 2, 3], both at weak [4, 5, 6, 7, 8] 
and strong coupling [9, 10, 11]. The computation of certain four-point functions at strong 
coupling has also been performed [12]-[25]. In this paper we use the integrability-based 
approach to correlation functions introduced in [7] to compute four-point functions of single 
trace gauge-invariant operators in the SU{2) sector of A/" = 4 SYM. Furthermore, motivated 
by the weak/strong coupling match for three-point functions in the classical limit reported in 
[8] , we will show that the same exact match occurs for the four-point function of two heavy 
non-BPS operators and two light BPS operators. In this introduction, we anticipate some 
of our main results and set our notations and conventions. 

Conformal invariance of A^ = 4 SYM fixes the two-point and three-point functions of its 
operators to take the form 

G2iXi,X2) = {0^ix,)di{x2)) 
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and 



A^|xi2|^i+^2-A3U23|A2+A3-Ai|3.g I 



where Xjj = Xj — Xj, Aj are the dimensions of the operators, C123 are the so-called structure 
constants and we have normalized the two-point functions to one. The structure constants are 
not completely unambiguous, since multiplying Oi by a phase will change C123 by that phase. 
On the other hand the absolute value IC123I is unambiguously defined. As opposed to the 
two- and three-point functions, the form of four-point functions is not uniquely determined 
in A/" = 4 SYM. Indeed, conformal invariance only tells us that the four-point function will 
depend on the cross-ratios 
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that we can form with the positions of the operators. Namely, in the planar limit 
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where A = Y,t=i ^i- 

At weak coupling, the tree-level computation of correlation functions is simply given by 
summing over all possible Wick contractions between the (constituents of the) single trace 
operators. Each single trace operator is an eigenvector of the one-loop dilatation operator. 
This is necessary to lift the degeneracy present at 't Hooft coupling A = 0, where there is a 
large number of operators with the same classical dimension [26]. The problem of performing 
the Wick contractions between these eigenvectors is purely combinatorial and can be tackled 
by exploting the integrability of the theory using the tools developed in [7] . 

We will consider single trace operators Oi made out of Lj scalars. Then, at tree level. 
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G2{xi,X2) = -. T^pr , G3(xi,a;2,X3) 
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where Uj is the number of Wick contractions between operators Oi and Oj. Note that given 
the lengths of the three operators the number of Wick contractions between the different 
single traces is uniquely fixed for three-point functions (2/12 = Li + L2 — I/3, etc). This is 
not the case for four-point functions. Instead^ 

G'4(Xi,X2,X3,X4) = ^ Y, rfllffe- ^^^ 

all possible {Uj} ^^ . '-^ 

As explained in this paper, each Ci234;{ii } can be computed using the integrability techniques 
developed in [7]. In other words, we can fix the function f{a,b) in (1) to leading order at 
weak coupling using integrability. 



"'^For non-planar diagrams, each Ci234;{/. } will come with additional 1/N factors. However, in this paper 
we will focus on planar diagrams, such that each Ci234;{/. } is simply a number. 
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Figure 1: Four-point function of SU{2) operators at tree level. All contractions are such that R- 
charge is preserved. The black (solid) lines represent vacuum, fields, while the red (dashed) lines 
represent excitations. Our setup is such that operators O3 and O4 do not interact among themselves. 
This restriction can be trivially relaxed at weak coupling using integrability techniques, see section 
5. 

We can consider four operators Oi such that only one choice of {Uj} survives due to 
R-charge conservation. This is the case for the setup considered in figure 1, for which our 
convention for the charges of the operators is the following 
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while the total lengths of the four operators will be denoted by 

Li = Ji + J2, L2 = Jl + J2-jl+J2-ki + k2, L-i=ji+J2 



h + k2. (4) 



For example I23 = J2, ^24 = ^2, etc. More importantly, l^^ = 0. We will demonstrate 
analytically and numerically that in the classical limit Li,L2 ^ L3,L4^ with Oi ~ O2, the 
four-point function factorizes as 
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which is somehow natural given figure 1. 

We will also consider more general four-point functions, for which different {lij} are 
allowed. In this case, we will explain that in the classical limit the term with I34 = in 
(2) dominates over the terms with l^^ 7^ 0.^ Furthermore for these more general four-point 
functions, we will show that in the classical limit 
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Cl234,{«,;,} ^ C123C124 , for /34 = . (6) 

The factorization (5) and (6) have interesting implications. A similar factorization of the four- 
point function of two large operators and two small operators was reported in [12] at strong 
coupling. In the Frolov-Tseytlin limit a match between weak and strong coupling for three- 
point functions was recently demonstrated in [8]. Hence, the factorization presented here 
extends this match to four-point functions of two large operators and two small operators. 
The numerical check alluded to above is extremely important in validating this match since 
we know that the issue of back-reaction is very important [8] . It could have been that further 
subtle issues arise for four-point functions. Our numerics indicate that this is not the case. 

The paper is structured as follows: in section 2 we review the holographic computation 
of the three- and four-point functions relevant to this paper. In section 3 we use integrability 
techniques to compute tree-level four-point functions at weak coupling. In section 4 we 
show numerically and analytically that in the Frolov-Tseytlin limit, the computations at 
weak and strong coupling of the four-point function of two heavy and two light operators 
match exactly. In section 5 we study a more general class of four-point function than the 
one considered in section 2 and study their classical limit using coherent states. Appendix 
A contains a short review of the main integrability tools to compute correlation functions 
in A/" = 4 SYM, while in appendix B we use integrability techniques to obtain a general 
formula for the four-point functions considered in section 5. Finally, appendix C contains 
the Mathematica codes needed to compute the main results of this paper. 



2 Strong coupling 

In [12], the holographic computation of four-point functions G'4(xi, X2, X3, X4) of single trace 
gauge-invariant operators was considered. It was argued that if the charges of the operators 
Oi and O2 are much larger than those of the operators O3 and O4, such that Oi ~ 025 then 
the four-point function factorizes into a product of two three-point functions 

CT4(^a;i, X2, X3, X/i) — — — r . {( } 

LT2ia^i5 X2) 

This is depicted in figure 2. If we take O^ and O4, to be light chiral primaries and consider 
the setup shown in (3), then only one choice of {Uj} survives and expression (7) simply 
translates into 

•-"1234 ~ '-"123'-"124 ) \° ) 



^We do predict the form of these other terms in the coherent state language. It would be very interesting 
to refine the strong coupling computation in order to probe these terms and verify if a match between strong 
and weak coupling holds for them as well. 




Figure 2: The strong coupling computation of the four-point function of two heavy and two light oper- 
ators factorizes into a product of two three-point functions, each of which involves a two-dimensional 
integration of a boundary-to-bulk propagator over the full string worldsheet. In the Frolov-Tseytlin 
limit both integrations become localized in the slice t^. = and only the integrals over a survive. In 
this way the computations at weak and strong coupling are matched. 

where • denotes a non-BPS operator, while o denotes a BPS operator. Consequently, we just 
need to know how to compute the holographic three-point function of two large operators 
and a light BPS operator. 

In order to compare the strong coupling results with the classical limit of our weak 
coupling computation, we need to take the Frolov-Tseytlin limit [27] of (8).'^ Luckily, we can 
simply use the results of [8] to compute the three-point functions appearing on the r.h.s of 
the factorization formula (8). Using the charges of the four operators as given in (3), at the 
end of the day one obtains (see section 2 of [8] for details) 
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■^Recall that given a classical string solution with total charge J, the Frolov-Tseytlin limit corresponds to 
taking A, J — ?► oo with A/J^ <C 1. For more details, see section 2 of [8]. 



C°24 is of course exactly the same as C^23 with ji -^ ki and Ui, U2 are determined by the 
classical string solution in S^ dual to operators Oi and 02- Recall that one takes the ratio 
of the structure constant of interest to the three-point function of three BPS operators in 
order to remove any dependence on the normalization convention. We also take the absolute 
value of the ratio to avoid the ambiguity that arises when multiplying each operator in the 
three-point function by a phase. Hence, the factorization formula (8) in the Frolov-Tseytlin 
limit reads 

n23'^124- (12) 



^"1234 

This is the main formula of this section 



'-'1234 



'-'123 '-'124 



A few remarks are in order. The charges of the operators in the four-point function in 
the factorization formula (7), or equivalently in (12), are those specified in (3). However, 
the charges of the operators entering the three-point functions in (7) and (12) must be, due 
to charge conservation, slightly different. This point was already raised in [12]; however, 
the operators considered therein were made out of the complex scalars Z and Z only. Since 
in this paper we are considering operators that also have X and X, let us give the correct 
prescription for the charges of the operators in the three-point functions entering the factor- 
ization formulas above. Given the charges shown in (3), the charges of the operators in the 
three-point function C123 are {^1,^2}, {-'^i — ji, -^2 + J2} and {ji, j^}, respectively. Similarly, 
the charges of the operators in C124 are { Ji, J2}, { ^i — A^i, J2 + ^2} and {fci, ^2}. 



3 Weak coupling 

In this section, we will describe the computation at weak coupling of tree-level four-point 
functions of generic operators in the SU{2) sector of A^ = 4 SYM. In the spin chain language, 
each single-trace operator Oi is represented by a Bethe state on a closed spin chain [28, 29] 
that we denote by |^I/i). For example, for an operator made out oi L — N complex scalars Z 
and A^ complex scalars X, we have 

l^)= Yl ^nu...,n^TT(^Z...ZXZ...ZXZ...y (13) 

where 4'ni,...,nM is the so-called Bethe wave function, whose specific form can be found in [7]. 
We can think of the scalars X as excitations propagating in the Z vacuum. These excitations 
are conveniently parametrized by their rapidities or Bethe roots Ui = ^ cot y, where Pi are 
the momenta of the excitations. 

We will consider the configuration in figure 3. Below we will show how to compute the 
quantity C1234 appearing in (2), where we dropped the subscript {lij} given that in our setup 
there is only one such choice, which is completely determined by the charges of the operators 
given again by (3). 

We should stress that in this section, we are not assuming that C3 and O4 are small op- 
erators. Hence, all the results from this section are valid for generic SU{2) operators obeying 
the setup shown in figure 3. Also, as we mentioned in the introduction, the configuration we 
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Figure 3: Setup for the computation of the four-point function at weak coupling. The black (solid) 
lines represent vacuum fields, while the red (dashed) lines represent excitations. The top figure 
shows the vacuum (Z or Z) and excitation choice (X or X) for each operator. The figure at the 
bottom shows the labelling of the excitations of each operator and the different partitions needed to 
perform the Wick contractions. 



are considering does not include interactions between operators O3 and O4. However, it is 
trivial to modify the formulas we present below to consider the more general case in which 
O3 and O4 interact with each other.'' 



3.1 Four-point functions by brute force 

We can always do a brute force computation to determine C1234 using the explicit form of 
the Bethe states representing each of the operators, see (13). Then, denoting the lengths of 
each operator by Li, see (4), we have 

Li— ji — fei min{i,Ji— ji — fei} 

cuu^n E E E E 

1=0 |ai|=0 l<ni<...<n|Qj^|<Z '^<rni<...<mj^_j^_f.^_^^^^^<Li—ji — ki—l 

X ?//^^ 

T^ni,...,n|„j^|,/+l,...,/+fci,fci+/+mi,...,/ci+/+mjj^_jj^_t.j^_|Q^|,Li-ji+l,...,Li 

X ^^^^ 

X 4%...,n ^lU^i,...,U (14) 

where ip^^"^ is the Bethe wave function associated to operator Oi, \ai\ is the number of Bethe 
roots in partition ai, see figure 3, and^ 



, L1L2L3L4 

This factor takes into account the equivalent ways of breaking the spin chains (i.e. due to 
cyclicity we can rotate each chain before cutting it). Finally, A/i denotes the norm of operator 
Oi, which is defined as 

^^= E (C..,n.J*(C..,n.J, 

l<7ii<---<njv^<Li 

where iVj is the number of excitations of operator Oi. 

Since each Bethe wave function ip^'^^ has N^l terms, we see that for large Ni, Li equation 
(14) is computationally extremely inefficient due to the huge number of terms involved. 
Below, we will see that using integrability techniques, we are able to simplify the computation 
of C1234 significantly. 



^We consider such four-point functions in section 5. 

■'^We can also include a symmetry factor in SI to take into account the case when O3 or O4 are dropped 
from the four-point function: 



fl — \ — - 3 ^ (15) 

V MAA2AA3AA4 (£l-Jl-fci)(e(i3)e(L4)-5L3>o5L,>o)+l' ^ ' 

with Q{x) being the Heaviside theta function, such that 8(a;) = 1, for x > and Q{x) = 0, for a; < 0. In 
this case we recover the brute force formula for three-point functions (see equation (85) of [7]). 
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3.2 Four-point functions from integrability 

The combinatorial problem associated with the multiple Wick contractions required to com- 
pute the tree-level four-point function can be solved using the integrability tools introduced 
in [7]. One simply needs to follow these steps: 

• We cut each of the closed spin chains 1 and 2 into four open subchains. Formally this 
means that we write each of the states |\l/i) and |\E'2) as a linear combination of tensor 
products of four states in open subchains. Schematically, i.e. leaving out the sums over 
the different partitions and the factors arising from cutting the chains (see appendix 
A), and using the notation in figure 3: 

l^i) -^ |ai) ® |tt2) ® las) ® las) , 1^2) ^ |/3i) ® |/32) ® IPs) ® IPs) ■ 

• Similarly, we cut each of the closed spin chains 3 and 4 into two open subchains, i.e. 
we write each of the states |\l/3) and |\1'4) as a linear combination of tensor products of 
two states in open subchains. Schematically 

1^3) -> |7) ® It) , 1^4) ^\6)®\6). 

• In order to perform the Wick contractions among the operators in the four-point func- 
tion, we first need to flip some subchain states. Again, leaving out the sum over par- 
titions and the relevant factors coming from the flipping procedure, we schematically 
have 

1^1) -^ \ai) ® 1^2) ® 1^3) ® ("si , 1^2) -^ |/3i) ® {i3;\ ® {i3;\ o {i3;\ . 
1^3) ^ I7)® (7*1, 1^4) ^\S)(^{6*\. 

• We contract, or sew, the different subchain states as shown in figure 3. This involves 
the computation of scalar products of Bethe states. To efficiently perform them, we use 
the new recursion relation for SU{2) scalar products derived in [7], which we quickly 
review in appendix A. 

• Finally, we normalize the states and sum over the distinct ways of breaking them. 

All the necessary tools to explicitly perform this computation are reviewed in appendix A, 
where all the notation used in the formula below is explained in detail. After all the dust 
has settled, we get 

^i —li — '^i -fctiai xaiai £0.20-2 £0202 

Ci234=fi Z2 ^ ^ ^ Yl ^Z-h-ki-l^kl ^TWW ~ 

1=0 0\VJax={u}02yja2=ax0zVJas=a2 p2VJI32={v} •'< 

l/32|=|a3l 
l/32|=|ai| 

fl32h f 1^2/32 fhh Az}{z} 

y as f03"3 {01303 £0303 {•"} 02 J J> J> {z} J> 

^ ^l+ji+l-l J> J< ^J2+l+l^L2-J2-l Av}{v} ^Li+l Az}{z} 

J< J< 

X {(3*,\a,){p;\a,){{z*}\a,){a;\{w}). (17) 
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In section 4, we will use formula (17) to compute numerically tree-level four-point func- 
tions that satisfy the setup of figure 3. We will then extrapolate our numerical results to the 
case in which the lengths of Oi and O2 go to infinity and see that in this particular limit, 
the weak coupling result factorizes as in (5) and hence matches the strong coupling result in 
the Frolov-Tseytlin limit (8) (or equivalently, (12)).^ 

Let us stress that the integrability-based formula (17) proves to be far more efficient than 
the brute force computation (14). For the reader's convenience, we present in appendix C 
the Mathematica codes needed to compute four-point functions using both formulas, as well 
as some specific examples showing how to use the codes. As noted in that appendix, even 
when the number of excitations and lengths of the operators are not too large, the brute 
force formula (14) becomes computationally much slower than formula (17). 



4 Weak/strong coupling match 

In this section we will show that in the classical limit, the four-point function C"34 of two 
large non-BPS operators and two small BPS operators computed at weak coupling matches 
the strong coupling result. First, we will provide some numerical evidence to support this 
claim. Namely, we will use the weak coupling result (17) to compute the four-point function 
of interest and then extrapolate numerically our results to the case when the lengths of the 
two large operators go to infinity. By doing so, we will obtain a numerical match with the 
strong coupling result evaluated using equation (12). We will then provide an analytic proof 
of this match by representing the two large operators by coherent states. 

4.1 Numerics: SU{2) folded string 

Let us consider the following operators in the four-point function represented in figure 3. 
Operators O^ and 04 are taken to be the following small BPS operator: 

0^ = 0^ = 2 Ti{ZZXX) + Tr{ZXZX) . (18) 

Operators Oi and O2 are taken to be large non-BPS operators dual to the folded string 
with unit mode number. This is very similar to the three-point function setup of two large 
operators and a light BPS operator considered in [8]. We can directly borrow the results 
from that paper to compute the structure constants appearing on the r.h.s of (12) in the 
Frolov-Tseytlin limit (we refer the reader to section 3 of [8] for details). We obtain 



^Note that a priori it is not obvious at all that (17) factorizes as in (8) when Li,L2 — > 00. However, in 
section 4 we provide an analytic proof of the weak/strong coupling match from another perspective, namely, 
by representing Oi and O2 by coherent states. 
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Figure J^: Two Bethe roots configurations for a folded string with a = 1/3 and Li = 30, 150, 
represented by the larger and smaller bullets respectively. The horizontal and vertical axes are the 
real and imaginary part of u/Li, where u is the rapidity of an excitation, related to its momentum 



by u 



cot I . The case Li = 30 corresponds to the maximum number of roots that we used to 



compute ri234 with formula (17). Even though this number of roots was not very large, we see that 
they lie nicely along the cuts formed by the clearly classical configuration Li = 150. 



where a is the filhng fraction of Oi and q is related to it by 

■Ji . E(g) 



a 



J1 + J2 



1 



K(g) 



(20) 



with E(g) = EllipticE[q] , K(g) = EllipticK[q] in Mathematica. Hence, given a filhng 
fraction a, we can find the corresponding value of q from (20) and plug everything into (19) 
to obtain a number, which will be the analytical prediction from strong coupling. The first 
column in table 1 has different such predictions for three different filling fractions. 

Having in mind the same setup considered in the previous paragraph, the goal of this 
subsection is to evaluate numerically our weak coupling formula (17) by keeping a fixed 
and increasing the length of this operator. In order to evaluate (17), all we need are the 
Bethe roots of the two large operators.^ Since these are dual to the folded string, their 
Bethe roots are distributed along two symmetric cuts, as shown in figure 4. In the setup we 
are considering, operator O2 has Ji — 4 Bethe roots, while Oi has Ji of them. It is then 
natural to ask where do we put the extra four roots of Oi. Following the discussion in [8], 
we choose these four roots to lie on the already existing two symmetric cuts defined by the 
Bethe roots of 02- Given our setup, this means that we add two roots to each of the cuts. 
As we are about to see, this is the weak coupling choice that reproduces the strong coupling 
computation. 

It is computationally easy to find the positions of a very large number of Bethe roots, see 
[30]. However, due to the sum over /, the many sums over partitions and the scalar products 



^For the light operators, recaU that to obtain a BPS state in the spin chain language, we should send all 
rapidities to infinity, which is equivalent to have all excitations with zero momentum. 
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appearing in it, equation (17) is quite non-trivial to evaluate numerically if we consider a 
large number of excitations.^ In practice we were only able to evaluate it for configurations 
with 4, 6, 8 and 10 Betlie roots in operator Oi using the Mathematica codes presented in 
appendix C. We collected this data in the following form {Li,ri234}, where ri234 is the ratio 
on the l.h.s of (12). For example, for filling fraction a = 1/3, we obtained 

data={{12 , . 669138111344} ,{18,0. 654008855014} , {24, . 639613251186} , 

{30,0.631580997967}}; 

We then fit this data in order to find the large Li asymptotics of ri234. Since for a fixed a 
we only had four points of data, the fit we used was of the form ri234 = ao + ai/Li + a2/ L\ + 
a^/L\. Given the data for a = 1/3 shown above, we can simply run the following code in 
Mathematica 

Fit [data, l/L"Range [0,3] ,L] 

to obtain oq = 0.627249. In the second column of table 1 we present the value of oq for 
three different filling fractions. We see that, with very small error, the numerical weak 
coupling results approach the analytical strong coupling predictions in the Frolov-Tseytlin 
limit. Figure 5 shows our list of data points, their fits and the analytical predictions. 



a 



Analytical prediction Numerical extrapolation Relative error 



0.639186 0.627249 1.868% 

0.734274 0.734578 0.042% 

0.789894 0.791947 0.259% 



Table 1: Numerical data obtained with the weak coupling formula (17) (evaluated with the Bethe 
roots of Oi and O2 lying on the same classical cuts) compared with the analytical prediction from 
strong coupling (19) for different values of the filling fraction of operator Oi. 

Recall from the discussion in [8] that the other two choices for the positions of the extra 
roots of Oi are i) to put them at finite positions outside the existing cuts of O2 and ii) to put 
them at infinity. This is the issue of hack-reaction addressed in [8] . The triangles in figure 5 
represent the data points we obtained by evaluating (17) with option i). As we can see, the 
extrapolation in this case does not match the strong coupling result. Hence, we confirmed 
numerically that the issue of back-reaction is also important for four-point functions.^ 

Before closing this subsection, let us make some comments regarding the accuracy of our 
numerical results. By looking at table 1 and figure 5, the skeptical reader might think that 
the results presented in this subsection are not conclusive enough to claim a weak/strong 
coupling match for four-point functions obeying the setup of figure 3. However, we should 



^Note that (17) is quite more involved than the three-point function formula derived in [7] and used in 
[8]. 

^Of course, we expect the case in which the extra roots are placed at infinity to also differ from the other 
two choices. 
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Figure 5: The bullets correspond to the numerical data obtained by evaluating the weak coupling 
formula (17) with the roots of Oi and O2 lying on the same classical cuts, for different values of 
a = 1/3, 1/4, 1/5, from bottom to top. In this case the fits (solid black curves) perfectly match the 
analytical prediction from strong coupling (dashed gray lines) obtained from (19). The triangles 
correspond to the data obtained evaluating (17) with a = 1/3 and by placing the extra roots of 
Oi at finite values outside the cuts of 02- Clearly, this case does not match the strong coupling 
prediction, confirming that the issue of back-reaction is also important for four-point functions. 

note again that we only used four points of data to perform each fit, compared to the seven 
points used for the numerical weak/strong coupling match of three-point functions of two 
large and one small operators [8]. In that case, the accuracy in the numerical match was in- 
deed much higher. However, had we only considered the first four points of data obtained in 
[8] to perform the fit for the three-point functions match, the relative error between the ana- 
lytical and numerical results for a = 1/3, 1/4, 1/5 would have been 2.122%, 0.612%, 0.381%, 
respectively. Hence, we see that the relative errors for the weak/strong coupling match of 
four-point functions shown in table 1 are in fact smaller than the relative errors in the match 
of three-point functions. We can be confident that adding more points of data to our fits 
would simply increase the accuracy of the numerical results even further, confirming that 
the weak/strong coupling match presented here is indeed exact. 



4.2 Analytics: four-point functions from coherent states 

Let us now prove analytically the match between weak and strong coupling for the four- 
point functions of interest. ^*^ We will use an alternative approach to the one presented in 
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We will follow closely the logic of section 3 of [8], to where we refer the reader for more details. 
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section 3 and used in the previous subsection. Namely, we will consider the operators Oi 
and O2 to be classical with their charges and their lengths taken to be much larger than the 
corresponding charges and lengths of O3 and O4,. In this limit, the exact Bethe states (13) 
are well approximated by coherent states [31] 



|Oi) = ■■■® u(^)^® u(^)^®... (21) 

{0,\ = ■■■®(u(^) ®(u(^) ®... (22) 



where £j = Li/2'K and at each site |u) = Ui\X) + U2\Z), such that (u|u) = ii • u = 1.^^ We 
note that these states do not depend on the Bethe roots and are completely disentangled. 
The operators O3 and O4 are chosen to be vacuum descendant states, corresponding to the 
BPS operators 

03 = M {Ti[X^'Z^''] + permutations) , (23) 

04 = M {Ti[X'''Z^^] + permutations) , (24) 

where the normalization constants A/3 and A^ are one over the square root of the number of 
distinct permutations 



Now consider the setup of figure 3. The operator C3 is Wick contracted with Oi and O2 
at sites^^ q,q — l, . . . . For a given q and /, the Wick contractions between the large operators 
Ci and O2 are given by 

l+q h2+q+ki 

^,,1 = n ^(^^2) ■ u(v/:i) n ^(^VA) ■ u(j7^i) • (27) 

i=q+l j=q+l+ki + l 

Let us consider the contractions between O3 or O4 with the classical operators. Given the 
present configuration, the only term of O^ which gives a nonzero contribution is Tr[Z-'^X-'^] 
and similarly for 04. For a fixed q and I, we denote the Wick contractions of O^ and O4 
with the classical operators by J'^i and JT./ , respectively. We have approximately 



Jlf^-^s<[- 



.3)....W|_U,U 



J-i,?-M4^(^ 



ri^).,^r.Mf ^ + ^ + ^ \^kJ q + l + l 



£1 / ' V '^ 



^^Note that to describe a given coherent state, we use the same notation Uj used in the strong coupUng 
discussion, see equations (9) and (10). This is because these quantities wiU be matched between weak and 
strong couphng. 

^^Due to the periodicity of the spin chain, the sites q = 0, —1, —2, . . . correspond to Li,Li — 1, Li — 2 . . . 
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Finally, we get 
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Li Li 



'-^1234 — 2-^ 2-^-^<i^^^q,l '-^q,l ' 
/=0 q=0 



{2t 



Following [8], we will now argue that Xq^i is one to leading order. For that we may use 
a gauge invariance of the absolute value of C*2^l: by multiplying each single site state by 
a phase \u^^^) — )■ e~'^'^^^^\u^^^), |C*234| must be invariant. We fix this gauge by choosing a 
conformal-like gauge, resembling the Virasoro constraint in the string theory side, which 
reads^'' u ■ d^u = 0. Then, with this gauge and using 



r _ r ki+ji-k2- J2 
L,2 — L,i-\ , 

Ztt 



(27) can be rewritten as 



(29) 



2ttI 

Lq^i ~ exp / — log ( u 



J2 + h - Ji - h 

u — — ad^u 

2,TiLo 



27r(/i2-0 



X exp 



da' 

2^ 



log u 



32 + k2-3i-ki ;^ 

lTlL2 



l + O 



c, 



With this simplification, the remaining part of the formula becomes 



'-^1234 — 



Ci 



u 



32 



c, 



/=0 q=0 
Jl!j2! 

(ji+j2-l)!V (^1 + ^2-1)! 



Ml 






u\- 



g + / + l 

Co 



k,\k2\ 



2tt 



2lT 



(30) 



L, ^< {a)uf (a) I L, ^u^^ {a')ul^ (a') 



where in the second line we have used (29) and the periodicity of the chain, which allowed 
us to factorize the double sum in the first line. Upon dividing this expression by C^23C*i: 
which can be read from (11), we exactly obtain the strong coupling result (12)! 



oo 
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5 More general four-point functions 

In the case we studied in the previous sections, operators O3 and C4 did not contract 
between themselves, see figure 3. Such a configuration was motivated by the fact that at 



^■^The upper limit in the sum over I should be Li — ji — ki. We replaced it by Li, since Li 3> ji, ^2- 
^^In the continuum limit, we introduce the variable a = 27r J and u*^"^ — >■ u(cr). 
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Figure 6: Although this is a tree-level diagram, it is suppressed by 1/Li with respect to the case of 
figure 2 since in this case only one integration over the worldsheet is performed. This is analogous 
to the weak coupling configuration of figure 7b. 



strong coupling these were taken to be light operators and their interactions are suppressed. 
Indeed, in the tree-level diagram of figure 6 only a single integration over the worldsheet is 
performed and hence it is suppressed by 1/Li with respect to the diagram of figure 2.^^. We 
can therefore address the question of whether or not we still have a match with the strong 
coupling result when we consider a weak coupling setup such as the one shown in figure 7, 
in which the light operators can interact. In particular, is the strong coupling suppression 
mentioned above manifest at weak coupling? 

For the setup we consider in this section we need to specify the number of excitations 
{X and X) and vacuum fields {Z and Z) contracted between operators O3 and O4. We will 
denote these by mi and 777,2, respectively. In this case the four-point function reads^^ 

mill { ji ,ki } min {J2 M } 
G4{Xi,X2,X3,X4) = —— 2_^ /_^ C'l234,{mi,m2} ^{mi,m2}(3^1? ^252:3, X4) (31) 



mi=0 



m,2=0 



where 



Q{mi,m2}['^li^2-, 3^3) "^i) 



1 



\xi2\'^'^^^~^^~^^~^'^^~^'^'^'^\xi2,\^'^^' 



2-"2)|Xi4p(^i- 



■mi) 



X 



(32) 



X 



\^{Ji-mi)\j-^A2{k2~m2)\j-^A2(mi+m2) 



F23 



We will now compute Gi{xi,X2,x^^X4) for the case where Ci and O2 are heavy operators 
and O^ and 04 are light chiral primary operators, using the coherent state language. In 



""^^By \/Li supressed we simply mean that the case of figure 2 is of order 0{L\), while that of figure 6 goes 
as 0{Li) precisely due to the number of integrations over the worldsheet involved in each case. 
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We recall that we are excluding the disconnected diagrams that may appear when {ji, J2} = {fci, ^2}- 
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Figure 7: More general four-point functions in the SU{2) sector, (a) It is clear that the diagrams 
with no contractions between O3 and 04^ are planar for any / > 0. (h) If there are contractions 
between O3 and O4, the only planar diagrams are those in which 1 = 0. When Oj, and O4, are taken 
to be much smaller than Oi and O2, this weak coupling configuration is analogous to the holographic 
four-point function depicted in figure 6. 



appendix B, equation (43), we provide the formula for G4{xi,X2,Xs,X4) obtained using the 
integrabihty tools, valid for any four operators obeying the setup of figure 7. 

In the classical limit, we again impose the conformal-like gauge such that the contractions 
between heavy operators are approximately one. The contributions from the contractions 
between the light operators O3 and O4 with the heavy operators are respectively given by 

ji3) ^ f^i + ^2\ ^ n-^, (q\ ,,-m, (q_ 

^(4) _ f^l + ^2\ ,r k,-m, ( Q + I + ^ \ -M-n., (Q + I + I 

^nl,m^ .mo — _ ■'^4 "l I ~ I "■ 



where q is the site of the heavy operators where Os is inserted and the binomial factor 
("^Irt"^^) ^^ ^^^ number of contractions between the light and heavy operators. 

For / 7^ and mi + 1712 7^ the diagrams are non-planar, and hence 1/A^ suppressed. 
Hence, in the planar limit, the four-point function is given by 



Giixi, X2, X3, X4) ~ -^^^ ^{0,0} (a;i, X2, X3, X4) Jg%^o Jq%fl + 



i=0 g=0 
min{ji,fei} min{J2,fc2} Li 

mi=0 012=0 g=0 

(33) 



where the first term takes into account the contribution of the diagrams for which {m-i, m,2} = 
{0, 0} and any I. The second term corresponds to the diagrams with / = and mi + 1112 7^ 0. 
Taking the continuum limit of this expression, we obtain 

2-K 2-K 

Gi{xi,X2,x^,Xi) ~ ^^{o,o}(a;i,a:2,a;3,a;4)A/'3A/'4 / Li— nf (o-)^^ (a) / Li—u\^{a')ul^{a') + 



min{ji,fci} mm{J2,fc2} / , \ 2 

1 ^^^ ^;^ r /-, / /mi+m2\ ., ., 

"^ Ar2 ^ ^ ^mi+m27^0t/{mi,m2}(a;i,a;2,a;3,X4) I I A/3A/4X 

mi=0 m2=0 ^ ^ '' 

2-K 



X / Li^^f-'"^(a)«i^-™^(a)«t^--(a)^2'^-"^^(a), 

ZTT 







(34) 



The second term is suppressed by 1/Li compared to the first one and thus, we recover the 
strong coupling result (12). In this way, we made manifest that the l/Li suppression due to 
the interaction of light modes in the string theory side, see figure 6, has a dual description 
in the gauge theory side in the classical limit. It would be very interesting to compare 
the second term in (34) with the result of a strong coupling computation of the tree-level 
four-point function of figure 6 and see if a match occurs for these results. 
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Finally, note that in order to evaluate (34), we need to specify the charges of the four 
operators according to (3), the values of mi, 1712 and the variables ui, U2. The latter depend 
on the classical string solution dual to the heavy operators Oi and O2 (e.g. see section 3 of 
[8] for their explicity expressions when the heavy operators are dual to the folded string). 



6 Conclusions and discussion 

This paper was devoted to the study of four-point functions of single trace gauge-invariant 
operators in the SU{2) sector of TV = 4 SYM. By representing each operator as an SU{2) 
spin chain state, we applied the integrability techniques introduced in [7] which allowed us 
to solve the problem to leading order at weak coupling for any four operators obeying the 
setup of figure 3. 

We then studied the classical limit of such correlation functions, focusing on the case 
in which the charges of two of the operators are taken to be much larger than those of the 
other two. We provided numerical evidence for the match of this weak coupling result with 
the strong coupling computation in the Frolov-Tseytlin limit. Furthermore, we managed to 
explain analytically this match, using spin chain coherent states to describe the two heavy 
operators. 

An important aspect for the numerical match was related to the difference between op- 
erators Oi and 02- Given that operators O^ and O4, were taken to be small, it is reasonable 
to expect that Oi should be approximately the complex conjugate of 02- In the classical 
limit, the Bethe roots of operator O2 organize into cuts. The slight difference between Oi 
and the complex conjugate of O2 comes from the position of the few extra roots of Oi with 
respect to those of 02- Given the cuts defined by the roots of O21 there are three possi- 
bilities for where to put the extra roots of Oi. they can be placed on the existing cuts, at 
finite position outside the existing cuts, or at infinity. Based on a similar analysis for the 
three-point function problem [8] , we chose the first option to obtain the referred weak/strong 
coupling match. We also showed that, just as in the case of three-point functions, the issue 
of back-reaction is important for four-point functions. That is, the weak coupling results 
obtained using equation (17) are sensitive to the different choices for the positions of the 
extra roots of Oi in such a way that only the choice we made for them matches with the 
strong coupling result (12). It would be very interesting to study the issue of back-reaction 
at strong coupling. 

We also considered a more general class of four-point functions in the SU{2) sector of 
the theory by allowing contractions among all four operators. When O3 and O^ are light 
operators, we found two types of contributions at weak coupling: the diagrams in figure 
7a give the leading contribution, while those in figure 7b are suppressed in the length of 
the heavy operators. Furthermore, by representing the two heavy operators by coherent 
states, we computed the exact form of these supressed terms. At strong coupling, there is an 
analogous behaviour. Two possible contributions at leading order are represented in figures 2 
and 6, with the former being the dominant one and the latter being suppressed in the length 
of the heavy operators since a single integration is performed over the worldsheet. It would 
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be very interesting to perform the strong coupling computation of the tree-level holographic 
four-point function depicted in figure 6 and see if a match between weak and strong coupling 
occurs in this case. 

Also, using integrability techniques, we showed how to compute the four-point function 
of figure 7 to leading order at weak coupling for operators of arbitrary size, see (43). We 
expect this formula and also (34) to capture part of the result for the one-loop four-point 
function at strong coupling. ^^ However, it would be important to first check this claim in 
the case of three-point functions, see the discussion in [7, 8]. 

We should stress that our weak coupling formulas (17) and (43) are valid for any four 
operators in the SU{2) sector . In particular, once an equivalent computation is performed at 
strong coupling^^, it would be interesting to make a comparison of both results for the cases 
when more than two operators in the four-point function are heavy and see if a weak/strong 
coupling match occurs. 

The strong coupling results presented in section 2 can be generalized to the SU{3) sector 
[8]. The same is not true at weak coupling, where the necessary integrability tools are not 
yet known, but such problem is currently under investigation [34]. Nevertheless, the match 
of the four-point functions of two large and two small operators persists in the SU{3) sector. 
This can be easily seen using SU{3) coherent states as in [8]. 

Finally, it was recently reported [35] that the match of tree-level structure constants in 
the classical limit is also valid in the SL{2) sector of the theory. It would be interesting to 
use the integrability tools for the SL{2) sector presented in appendix A of [7] to see whether 
or not the match for four-point functions and the issue of back-reaction also extend to that 
sector. 
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^^As pointed out in [7], we need to do two things in order to get one-loop data from the tree-level formulas 
obtained using integrability. Firstly, we need to consider the 0{\) correction to the Bethe wave function of the 
operators, obtained by diagonalizing the two-loop Hamiltonian. Secondly, we have to consider Hamiltonian 
insertions [4, 5, 6] when performing the contractions between all operators. 

^^Recently there has been progress towards tackling the computation of general holographic three-point 
functions at strong coupling [8, 32, 33]. 
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A Cutting, flipping and sewing spin chains 

In this appendix, we quickly summarize the integrabihty tools introduced in [7] for the 
computation of correlation functions of single trace gauge-invariant operators in the SU{2) 
sector of A/" = 4 SYM. Note that all formulas we present in this appendix are in the coordinate 
Bethe ansatz base (see [7] for details). 

Cutting 

Consider a closed SU{2) spin chain described by a Bethe state parametrized by a set of A^ 
Bethe roots {u}. We can break it into a left and a right open subchain of lengths / and r as 



eaa 



IW)= E ^f WW \<^)^®\<^)r^ (35) 

ayja={u] J < 

where the sum runs over all 2^ possible ways of splitting the rapidities into two partitions 
a and a. For example, if A^ = 2, the possible partitions (a, a) would be ({}, {'Ui,M2}) , 
({■^i}) {""2}) , ({^2}, {wi}) , ({"Ui, U2}, {}). We are also using the following shorthand notation 






and 



Vj £ a i < j 



f{u) = 1 + - . (37) 



Flipping 

Consider again a Bethe state on a spin chain of length L. We can flip the state as follows: 

^°IW) = 4?iC^(K}|C', (38) 

where C stands for a charge conjugation which exchanges Z -H- Z and X ^ X in the operator 
language. We are using a shorthand notation in the same spirit as equations (36), where we 
would simply need to replace a by {u}. 

Sewing 

Consider two Bethe states on spin chains of length L parametrized by their Bethe roots Ui 
and Vi, with i = 1, . . . ,N. The scalar product S'Ar({f }, {u}) = {{v*}\{u}) can be efficiently 
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computed using the new recursion relation for SU(2) scalar products derived in [7]. It reads 

Sn {{Vl, . . . , Vn}, {Ui, ..., Un}) = '^K Sn-1 ({f 1, ...,Vn,---, Vn}, {Ul, ^2, • • • , Un}) 
- ^ Cn,m Sn-1 {{Ul, Vi, . . . Vn, ■ ■ ■ , Vm, ■ ■ ■ Vn}, {«!, ^2, • • • , Un}) AtoC ({«}, {f }) , (39) 

where a hatted Bethe root means that it is omitted and 
AtoC {{u}, {v}) ^ 



^{u}^{v}g{u+l}g{v-^^}f{n}{n} fMM ' 

N 
bn = g{Ul - Vn)a{Vn)d{ui) J]^ f {Ui - Vj)f{Vj - Vn) + («! ^ Vn) , 

N 
Cn,m = g{Ul - Vn) g{Ui - Vm) a{Vm)d{Vn) f{Vn - Vm) JJ /(w„ - Vj) f{Vj - Vm) + {u ^ Ul) , 

where 

g{u) = '- , a{u) = (^ + ' d{u) = [^ ' '^ ■ W 

Equation (39) provides a complete solution for any SU(2) scalar product and is the one 
we used to perform the explicit numerical evaluation of the four-point function at weak 
coupling, see equation (17). However, when one of the states in the scalar product is a 
BPS state, the formula for the scalar product simplifies dramatically [7]. Recall that in the 
spin chain language a BPS state is that in which all excitations have zero momentum or, 
equivalently, infinity rapidity. Hence, if we send the Bethe roots {v} to infinity in the scalar 
product above, we obtain the following simpler formula: 

<{->iw> = iTi™ 5: {-i^'ein. (41) 

y J< a.\Ja={u} 

B Four-point functions of section 5 using integrability 

We compute G^i^Xi, X2, x^, x^) for the setup of figure 7 using the integrability tools of appendix 
A. The steps we need to follow are very similar to the ones outlined for the setup of figure 
3, see section 3, except that now we also have to cut operators O3 and 04 more than once 
and perform some extra scalar products. In the planar limit, G^i^Xi, X2, X3, X4) is given by 

_ Li-Ls+ji-fci 

G4,{Xi,X2,X3,Xi) =—Q{ofi}{Xi,X2,X3,X4) ^ C'l234,{0,0,i} + 

1=0 

min{jr'i,fci} min{J2,fc2} 
~'~ 1V2 Z-^ Z-^ Om^+ma^O ^{mi,m2}(2^1)2;2,3^3)2;4) W234,{mi,m2,0} ) 



rrai=0 r?i2=0 



(42) 
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Figure 8: Setup for the computation of the four-point functions of section 5 at weak coupling. The 
black (solid) lines represent vacuum fields, while the red (dashed) lines represent excitations. The 
figure shows the labelling of the excitations of each operator and the different partitions needed to 
perform the Wick contractions. 



where 

Cl23i,{mi,m2,l} 
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aiUai={u} 02Ua2=ai (5iUl3i={v} /32U/32=/3i 72U72={ty} S2US2={z} 

IfeNlasI l72|=|/3i| |52|=|a2| 
l/32|=|«l| 



■,"1 



X F € 

Li-L3+ji-ki+mi+m2~l fci-mi 



foioi faiai fa2a2 fci2<^2 fci2Ci2 
(52 J J< J J< J< 



Au}{u} 



.^1 J2 J2 



fh(5i iPiPi rl32l32 f 132(^2 rl32l32 



j{v}{v} 



^L3-ji-m2 ^h+m2+l 



fy2f2 fl2"l2 ^7272 
Aw}{w} 



p{^} P&2 J2 

L4— fci— m2 mi+m2+l fei+m2+l 



^<52<52 ^^2^2 i'<52(52 
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(43) 
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where Q is given in (16) and the labehng for the Bethe roots is indicated in figure 8.^^ The 
first term in (42) corresponds to the contribution of diagrams for which mi,m2 = and 
/ > 0. In this case, expression (43) simphfies since the partitions 62 and 72 become empty. 
The second term in (42) takes into account the contribution of planar diagrams where the 
hght operators can contract and therefore we should have I = 0. Equation (43) is valid for 
any four operators obeying the setup presented in figure 7. 



C Mathematica codes 

In this appendix we provide the Mathematica codes for computing the four-point func- 
tions obtained in section 3. We implement both the brute force result (14) as well as the 
integrability-based formula (17) and give some explicit examples at the end of this appendix. 
With some slight modifications of the code presented below, equation (43) can be imple- 
mented in a similar manner, but we leave that to the interested reader. 

Some useful functions 

Le=Length; 

f[u_]=l+I/u; g[u_]=I/u; n[0]=0; m[0]=0; 

f [ll_List , 12_List] : =Product [f [11 [ [ j 1] ] -12 [ [ j 2] ] ] , { j 1 , Le [11] } , { j 2 , Le [12] }] 

h [ll.List , 12_List] : =Product [h [11 [ [ j 1] ] -12 [ [ j 2] ] ] , { j 1 , Le [11] } , { j 2 , Le [12] }] 

fs[ll_List]:=Product[f [ll[[jl]]-ll[[j2]]],-[jl,Le[ll]},{j2,jl+l,Le[ll]}] 

gs[ll_List]:=Product[g[ll[[jl]]-ll[[j2]]],{jl,Le[ll]>,{j2,jl+l,Le[ll]}] 

fb[ll_List]:=Product[f [ll[[jl]]-ll[[j2]]],{jl,Le[ll]>,{j2,jl-l>] 

gb[ll_List]:=Product[g[ll[[jl]]-ll[[j2]]],{jl,Le[ll]>,{j2,jl-l>] 

gp[l_List] :=Times@@g [1+1/2] ; gin[l_List] :=Times@@g [1-1/2] 

e [l_List] : =Times@@( (1+1/2) / (1-1/2) ) 

Dvd [ls_List] : = ({Complement [Is , #] , #>&) /©Subsets [Is , {0 , Le [Is] }] ; 

BPS[N_] :=Table[10-(10+j),{j,N>] ; 

Four-point functions by brute force (14) 

Off [Det: :matsq] ; 

S[x_,y_] :=(x-y+I)/(x-y-I) 

Wave [l_List] : =Block [{p=Permutations [Range [Le [1] ] ] , i , j ]■ , Sum [A [p [ [i] ] ] 

Product [((l[[p[[i,j]]]]+I/2)/(l[[p[[i,j]]]] -1/2)) -n[j],{j,l,Le[l]}],{i,l,Le[p]>] 

//.{A[{a___,b_,c_,d___}] :>S[l[[b]] ,l[[c]]]A[{a,c,b,d>]/;b>c} 



^^Just like for the setup of section 3, we can also include a symmetry factor in XI to take into account the 
case when O3 or O4 are dropped from figure 7. In this case, it reads 



-S 



Lii2i3 ^^"^4 ^^" 1 



^irii ,m2 



AA1AA2AA3AA4 (Li - L3+J1 - fci + mi + m2)(e(L3)e(L4) - Sl,>oSl,>o) + 1 
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/ . {A [a ] : >1/ ; a==Range [Le [a] ] >] ; 

dphi [L_ , l_List] : =Det@Table [-If [i== j , L/ (1 [ [i] ] '2+1/4) -Sum [2/ ( (1 [ [i] ] -1 [ [k] ] ) "2+1) , 

{k, Le [1] >], 0] -2/(l+(l [ [i]]-l[[j]]) -2), {i,Le[l]},{j,Le[l]}]/.Det [{}]->! 

normdet [L_ , l_List] : =Block [{inf , f in, onlyf in} , inf =Select [1 , Abs [#] >10''8&] ; 

If [inf =={} , f in=l , f in=Select [1 , MemberQ [inf , #] ==False&] ] ; 

onlyf in=(fb [fin] dphi [L,f in] )/(fb[Conjugate[f in]]gp[f in]gm[f in] ) ; 

If [inf=={},onlyfin, (((L-2Le[fin]) !Le[inf] !) onlyf in)/ (L-2Le [fin] -Le [inf] ) !]] 

ClearAll [C1234bf ] 

C1234bf [Ll_ , Nl_ , L2_ , N2_ , L3_ , N3_ , L4_ , N4_ , ll_List , 12_List , 13_List , 14_List] : = 

Block [{j ,k,psis, norms, lini,lim2},psis [r_,s_] :=(Wave [11] /.{n[j_] :>r+j-s/;s<j<=s+N4, 

n [k_] : >N4+r+m [k-s-N4] / ; s+N4<k<=Nl-N3 , n [q_] : >Ll+q-Nl/ ; q>Nl-N3}) 

(Wave [12] / . {n [k_] : >L2- (L4-N4+r+m [Nl-N3-N4-s-k+l] ) +1/ ; k<=Nl- (N3+N4+s) , 

n [j _] : >L2-n [Nl- (N3+N4) - j+1] +1/ ; j >N1- (s+N3+N4) » (Wave [13] / . n [ j _] -> j ) 

(Wave [14] / . n [ j _] ->L4-N4+ j ) ; 

norms=normdet [LI , 11] normdet [L2 , 12] normdet [L3 , 13] normdet [L4 , 14] ; 

lim[s_,r_] :=Sequence@@Table[{n[j] ,n[j-l]+l,r},{j ,s>] ; 

lim2[s_,r_] :=Sequence@@Table [{m[k] ,m[k-l]+l,Ll-N3-N4-r},{k,Nl-N3-N4-s}] ; 

Sqrt[Ll L2 If [L3==0, 1 ,L3] If [L4==0, 1 ,L4] /norms] 

Sum[Sum[If [L3==0 | | L4==0, 1/(L1-N3-M4+1) , l]psis [r,s] ,lim2[s,r] ,lim[s,r]] , 

{r , , Ll-N3-N4> , {s , , Min [r , N1-N3-N4] >] ] 



Four-point functions from integrability (17) 

AtoC[u_List,v_List,L_] :=Block[{a,d,f ,g,z>,a[z_]=(z + I/2)~L;d[z_]=(z-I/2)~L; 
f[z_]=(z + I)/z;g[z_]=I/z; (l/(Times@@g[u+I/2]Times@@g[v-I/2] 
TimesOSa [v] Times@@d [u] ) ) ( 1/ (Product [f [u [ [ j 1] ] -u [ [ j 2] ] ] , -[ j 1 , Le [u] } , 
{j2,jl+l,Le[u]}]Product[f [v[[jl]]-v[[j2]]],{jl,Le[u]>,{j2,jl-l>]))] 
ClearAll [NewrecAl] ; 

NewrecAl [u_ , v_ , L_] : =NewrecAl [u , v , L] =Module [{N=Le [u] ,n,i,j,a,d,z,f,g,nl, n2} , 
a[z_] = (z+I/2)-L; d[z_] = (z-I/2)-L; f [z_]=(z+I)/z; g[z_]=I/z; 
Suin[g[u[[l]]-v[[n]]](Product[If [j!=n,f [u[[l]]-v[[j]]],l],{j,l,N>] 
Product[If [j!=n,f [v[[j]]-v[[n]]],l],{j,l,N}]a[v[[n]]]d[u[[l]]]- 
Product[If [j!=n,f [v[[j]]-u[[l]]] ,1] ,{j,l,N}] 

Product [If [j!=n,f[v[[n]]-v[[j]]],l],{j,l,N}]a[u[[l]]]d[v[[n]]]) 
NewrecAl [Drop [u , 1] , Drop [v , {n}] , L] , {n , 1 , N>] - 
If [Le[v]==l,0,Sum[(g[u[[l]]-v[[nl]]]g[u[[l]]-v[[n2]]]) 
(Product [f [v [ [nl] ] - j ] f [ j -v [ [n2] ] ] , { j , Join [Take [v , nl-1] , 
Take[v,{nl+l,n2-l}] ,Take[v,{n2+l,Le[v]}]]>]f [v[[nl]] - v[[n2]]] 
a[v[[n2]]]d[v[[nl]]]+Product[f [j-v[[nl]]]f [v[[n2]]-j], 
{j, Join [Take [v, nl-1] ,Take[v,{nl+l,n2-l}] , 

Take[v,{n2+l,Le[v]>]]}]f[v[[n2]]-v[[nl]]]a[v[[nl]]] d[v[[n2]]]) 
NewrecAl [Drop [u,l] , Join [{u[[l]] }, Join [Take [v, nl-1] ,Take [v,{nl+l,n2-l>] , 
Take[v,{n2+l,Le[v]>]]] ,L] ,{n2,2,N},{nl, l,n2-l>]]] 
NewrecAl [{}, {}, L_] = 1; 
SProduct [u_List , v_List , L_] : =AtoC [u , v , L] NewrecAl [u , v , L] 
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ClearAll[C1234int] 

C1234int[le_,Ll_,Nl_,L2_,N2_,L3_,N3_,L4_,N4_,ll-List,12_List,13_List,14_List] := 

Block [{al , alb , a2 , a2b , a3 , a3b , bl , bib , b2 , b2b , b3 , b3b , dvl , dv2 , dv3 , dv4 , norms} , 

dvl=Dvd[ll]; dv2=Dvd[12] ; 

norms=nornidet [LI , ll]normdet [L2, 12]iiormdet [L3,13]iiormdet [L4,14] ; 

Sqrt[(Ll L2 If [L3==0, 1 ,L3] If [L4==0, 1 ,L4] )/norms] 

SumElf [L4==0| |L3==0,1/(L1-N3-N4+1),1] (al=dvl [ [i , 1] ] ;alb=dvl [[1,2]] ; 

b2=dv2 [ [ j , 1] ] ; b2b=dv2 [ [ j , 2] ] ; If [Le [al] ==Le [b2b] , dv3=Dvd [alb] ; 

Sum [a2=dv3 [ [k , 1] ] ; a2b=dv3 [ [k , 2] ] ; If [Le [a2] ==N4 , dv4=Dvd [a2b] ; Sum [a3=dv4 [[1,1]]; 

a3b=dv4 [ [1 , 2] ] ; If [Le [a3b] ==N3&&Le [a3] ==Le [b2] , (1/ (f s [11] fs [12] ) ) 

e [alb] - (Ll-N3-N4-le) e [a2b] "N4 e [a3b] " (N3+le+l)e [12] " (L3-N3+le+l) 

e [b2b] - (L2-L3+N3-le) e [14] " (L4+1) f [al , alb] f s [al] f [a2 , a2b] f s [a2] f [a3 , a3b] f b [a3b] 

f s [a3] f [b2 , b2b] f b [b2] f b [b2b] (f b [14] /f s [14] ) SProduct [a3 , b2 , le] 

SProduct [al , b2b , Ll-N3-N4-le] SProduct [a2 , 14 , N4] SProduct [13 , a3b , N3] , 0] , 

{l,l,Le[dv4]}],0],{k,l,Le[dv3]>],0]),{i,l,Le[dvl]},{j,l,Le[dv2]>]]; 

C1234int [Ll_ , Nl_ , L2_ , N2_ , L3_ , N3_ , L4_ , N4_ , ll_List , 12_List , 13_List , 14_List] : = 

(temp[le_]=C1234int[le,Ll,Nl,L2,N2,L3,N3,L4,N4,ll,12,13,14] ; Sum [temp [le] , 

{le,0,Ll-N3-N4}]) 



Examples 

Let us now show how to use the code to compute four-point functions. Note that the functions 
C1234bf and C1234int can be used to compute any four-point function obeying the setup 
of figure 3. In the examples below, us, vs, ws and zs are the Bethe roots of operators Oi, 
O2, O3 and O4, respectively, for the charges indicated in each case.^° Let us consider first 
an example in which all four operators are non-BPS. Running the following code 

Ll=12 ; Nl=6 ; L2=12 ; N2=2 ; L3=4 ; N3=2 ; L4=4 ; N4=2 ; 

us=-[0. 676245041405- 

. 993633391204 1,0. 678017442247 , . 676245041405+0 . 993633391204 I , 

-0 . 676245041405+0 . 993633391204 1,-0. 678017442247 , -0 . 676245041405- 

0.993633391204 I}; 

vs=-[l . 702843619444 , -1 . 702843619444} ; 

ws={0 . 288675134594 , -0 . 288675134594} ; 

zs={0 . 288675134594 , -0 . 288675134594} ; 

C1234bf [Ll,Nl,L2,N2,L3,N3,L4,N4,us,vs,ws,zs] 

C1234iiit [LI , Nl , L2 , N2 , L3 , N3 , L4 , N4 , us , vs , ws , zs] 

we obtain a perfect agreement between the brute force computation and the integrability- 
based formula, giving 

^•••; = 1.28031 - 0.66373 i . 

As we mentioned in the main text of the paper, the integrability-based formula is computa- 
tionally much faster than the brute force formula. This can be checked in Mathematica for 



^'^ Again, it is computationally easy to find the positions of a very large number of Bethe roots, see section 
7 of [30]. 
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the example above by simply using the function AbsoluteTiming when executing C1234bf 
and C1234int. The result is that the former takes 18 seconds to compute, while the latter 
only takes 1 second! Of course, the gain in efficiency is even more notorious as we increase 
the length and number of excitations of each operator in the four-point function. 

Now, let us consider a case in which O3 and O^ are BPS operators. The following 
configuration corresponds to the second point of data for a = 1/3 used in section 4. Running 

Ll=18 ; Nl=6 ; L2=18 ; N2=2 ; L3=4 ; N3=2 ; L4=4 ; N4=2 ; 

us={l. 892118528217- 

1 . 114121041954 1,1. 976939483361 , 1 . 892118528217+1 . 114121041954 I , 

-1 . 892118528217+1 . 114121041954 1,-1. 976939483361 , -1 . 892118528217- 

1.114121041954 I}; 

vs={2 . 674763752754 , -2 . 674763752754} ; 

ws={10-10, 10-15}; 

zs={10-ll, 10-16}; 

C1234int [LI , Nl , L2 , N2 , L3 , N3 , L4 , N4 , us , vs , ws , zs] 



we obtain 



C''°l = -5.38224 - 6.23581 i . 



Upon diving this result by the product 0^23(^124,^ which can be implemented using the code 
presented in appendix D of [7], and taking the absolute value of the resulting ratio, we get 
^"1234 = 0.65401, which is exactly the second blue bullet plotted in figure 5. 

Finally, we should mention that when some of the operators in the four-point function 
are BPS, it is computationally more efficient to replace the appropriate general scalar prod- 
ucts SProduct (which compute formula (39)) appearing in C1234int by a function that 
implements the scalar product between a Bethe state and a BPS state (41). 
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